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(54) METHOD AND APPARATUS FOR ELECTROCHEMICAL MEASUREMENT USING 
STATISTICAL TECHNIQUE 


(57) A predetermined voltage is applied to a biosen- 
sor twice to promote an electrochemical reaction, and 
the following parameters (P1 and P2) are calculated 
from the values of detected current. A statistical tech- 
nique is used with these parameters to compensate for 
errors so that the concentration of an object can be 
determined. P1 : the ratio (If/lb) of the maximum current 
(If) or a current occurring the maximum in the first exci- 
tation to a current (lb) read at any point in the second 
excitation. P2: a current (lb) read at any point in the sec- 
ond excitation. 
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[0001] This invention relates to a method for electrochemically measuring a concentration of an analyte in a sam- 
ple. More specifically, this invention relates to a method for electrochemical measurement to obtain a concentration of 
analyte in a sample by correcting errors caused by sample physical properties and sensor sensitivity. 

BACKGROUND ART 

[0002] Recently in the fields of biochemistry, clinical medicine or the like, electrochemical measurement has been 
used to measure a concentration of an analyte in a sample in a rapid and simple manner. An electrochemical measure- 
ment is a method to measure an analyte in a sample by using an electric signal such as a current that is obtained from 
a chemical reaction or an enzymatic reaction. 

[0003] For example, the concentration of glucose in a blood is measured in the following process. Glucose as an 
analyte in a blood sample reacts with glucoseoxidase (GOD) to generate a GOD • H2 complex, from which electrons 
are liberated by an electron carrier such as potassium ferricyanide. The free electrons are captured in an electrode to 
calculate a current value, and thus, a glucose concentration is obtained by calculation from a calibration curve a cali- 
bration formula or the like that is produced previously. 

[0004] In an electrochemical measurement, a disposable device called a biosensor is commonly used. A biosensor 
has a reaction layer and an electrode system. 

[0005] A reaction layer comprises a reaction reagent, an enzyme or a matrix to react specifically with an analyte in 
a sample, and provides a reaction site. 

[0006] An electrode system comprises a working electrode and a counter electrode, and by applying voltage to per- 
form an oxidation-reduction reaction, the electrode system captures electrons that are generated by the chemical reac- 
tion occurring in the reaction layer, where the electrons are electric signals of a current from an electron carrier to the 
electrode. 

[0007] A biosensor is used in a combination with a measurement apparatus having various functions such as pro- 
viding a certain voltage in a predetermined period, measuring electric signals transferred from the biosensor, and con- 
verting the electric signals into a concentration of the analyte. Such a system is called a biosensor system. 
[0008] Among various methods to apply voltage to an electrode system, a method to apply voltage like a rectangu- 
lar wave with respect to time is called a potential step method. In a typical biosensor system, sample feeding to the bio- 
sensor is detected and subsequently, an open circuit or a voltage with substantially no current flow is provided to 
promote a chemical reaction. Alter a predetermined period of time, a fixed voltage is applied to deliver electrons 
between the electron carrier and the electrode, i.e., to perform an oxidation-reduction reaction. A state where the bio- 
sensor is provided with a voltage to perform a desired oxidation-reduction reaction is called an excitation state. In gen- 
eral, a current value at an arbitrarily-predetermined point of time during the excitation state is measured, and the current 
value is converted into a concentration of the analyte, using a calibration curve and a calibration formula that are pro- 
duced previously. 

[0009] Methods in which plural excitations and plural current measurements are performed are disclosed for exam- 
ple, in JP-2651278, JP-A-8-304340, and JP-A-1 0-1 01 30. 

[001 0] JP-2651 278 provides a method to determine whether a current flowing in a reaction site changes in accord- 
ance with a relationship with a certain Cottrell current 

[001 1 ] JP-A-8-304340 suggests reduction of measurement errors caused by a reduction type intermediate product 
that is generated during a storage of the sensor. 

[0012] JP-A-1 0-1 01 30 provides a method for discriminating a whole blood sample and an aqueous solution as a 
control in a biosensor system. 

[001 3] Measurement results obtainable by using the biosensor can include errors due to various factors. One of the 
factors causing such errors is a sample physical property. 

[0014] For example, when the sample is whole blood and the analyte is glucose in the blood, the hematocrit value 
(Hct) as a volume ratio of erythrocyte to the whole blood is found to cause errors in the measurement result, and differ- 
ences between individuals are great. The reason appears to be a rising viscosity of the sample. 
[001 5] Concentrations of neutral fat and protein in the blood also affect the measurement result. In a measurement 
of an analyte in a blood, sample physical properties, such as blood cells, lipid and protein, will cause measurement 
errors. 

[0016] Generally in a conventional method to avoid any influences by such errors, the composition of a chemical 
reaction layer or of an electrode system in a biosensor is improved. For example, JP-A-62-64940 discloses a biosensor 
in which an enzyme to detoxify a measurement-interfering substance is immobilized on an enzyme-immobilizing mem- 
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brane. JP-A-61-304S discloses a biosensor comprising not only a biocatatytic electrode but an electrode to detect the 
quantity of measurement-interfering substances. JP-A-60-21 1350 discloses a biosensor comprising two electrode sys- 
tems including an electrode system containing an enzyme and a pigment and also an electrode system provided with 
a porous material containing a pigment only. However, these biosensors with complicated structures require compli- 
5 cated manufacturing processes, and the manufacturing cost is also raised. 

[001 7] Sensor sensitivity can be another factor causing measurement errors. In many cases, sensor sensitivity var- 
ies from one manufacturing lot of biosensors to another. Manufactures make sensor sensitivity correction tips for the 
respective lots, and ship the tips with their biosensors. 

[0018] When a lot for a biosensor is changed, a user should correct sensitivity by using a correction tip correspond- 
10 ing to the changed lot before he performs an ordinary measurement (JP-A-4-357452). 

[001 9] However, such an operation will impose extra work on the user. Moreover, since the range corrected with 
such a correction tip is determined based on the sensitivity during manufacture of the sensor, changes of the sensor 
sensitivity over time, which occur after shipping, will not be corrected. 

[0020] Since an electrochemical measurement includes a chemical reaction, it is affected also by an environmental 
is temperature and a sample temperature. Namely, an environmental temperature and sample temperature also can be 
factors of measurement errors. JP-2748196 provides a method of correcting nonlinear temperature dependency of a 
chemical sensor. 

[0021] Therefore, the purpose of the present invention is to provide an electrochemical measurement method that 
can decrease work for users and manufacturers without requiring any complicated structure or process of manufactur- 
20 ing a biosensor and a measurement apparatus. Such a method can provide highly accurate results by correcting errors 
of concentration of an analyte in a sample. 

DISCLOSURE OF THE INVENTION 

25 [0022] To achieve the purpose, the present invention provides an electrochemical measurement to measure a con- 
centration of an analyte in a sample by using a biosensor having an electrode system and a chemical reaction layer. 
The measurement comprises calculation to obtain as parameters a current value obtained as a result of application of 
a fixed voltage to the sensor after feeding the sample and the ratio of the current value, and calculation of the concen- 
tration of the analyte by using a statistical technique. The parameters are set adequately for each error factor affecting 

30 the measurement results, and error factors can be corrected directly by selecting parameters that will be affected 
greatly by the sample physical properties, or by selecting current values that will be affected greatly by the sensor sen- 
sitivity. 

[0023] For the above-mentioned parameters, for example, the following parameters Pt and P2 are preferably used, 
and the parameters are obtained by applying a predetermined voltage twice to a biosensor after feeding a sample in 
35 order to promote an electrochemical reaction. 

P1 : a ratio (If/lb), where (If) is the value of maximum current or a current occurring after the maximum in a first exci- 
tation, and (lb) is the value of current (lb) read at any point in a second excitation 
P2: a current value (lb) read at any point in the second excitation 
40 With a statistical technique using these two parameters errors caused by sample physical properties and sen- 

sor sensitivity can be conected with high accuracy. Therefore, the present invention provides measurement with 
high reliability without requiring any complicated structures or methods of manufacturing sensors, or imposing extra 
work on users. 

In addition to the above-mentioned parameters Pt and P2, the following parameters P3 and P4 also can be 
45 used. 

P3: a value (l/AI(y)) obtained by normalizing a derivative value or a differential value of a current at any point in a 
second excitation, with a current value at the same point. 

Here, "normalization" means adaptation of a ratio of either the derivative value or the differential value to a cur- 
rent value at the same point of time in order to make the derivative value or the differential value a parameter not 
so dependent on the concentration of the analyte. When a ratio of the differential value to a current value is adapted, 
the current value can be read at any point between two points having a finite difference. 
P4: a ratio (lb(a)/lb(P)) of an initial current value (lb(a)) to a terminal current value (lb(P)) in a second excitation 

[0024] The inventors consider that these parameters are indices as shown below. However, as it still remains in the 
55 realm of speculation, one parameter will not always correspond to one index. 

P1 : mainly, a current value affected greatly by sample physical properties 

P2: mainly, a current value to indicate a concentration of an analyte in a sample 
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x rs So pS ns di,,usion of a compound and a mMute ° f a samp,e and a ~« * 

P4: mainly, a current value affected greatly by the electrode sensitivity of a sensor 

5 [0025] In the present invention, the current value lb at any point in the second excitation, which is used for calculat 
« Pa T ,erS "I "I" 2 .' iS Preferably me te ™ na " ~ -» W» - the second excitalon 

re analyuT 9 P ' nVent, ° n * 3 Parame,er ran9e that iS expec,ed from 1,16 concentration of an 

[0027] It is preferable in the present invention that correction formulas including a plurality of statistical techniaues 

M^ST^ IT" ° f 3 SamP ' e are Pr6pared and a " optimum ctSn 0 ' 
M^r^JZZZ m T 0 TH Perf0rm a " adeqUa,e C ° rreC,i0n ^^nding to either the environmental tem- 
perature or the sample temperature. The reason is that the current value obtained by applying voltage to the sensor to 

Var,6S COnSiderab ' y ^ ^ em '"° M ,empera,Ure °< >™Z»* temperature! and the 
ESL,^^ 8 H. the environmental '^perature or the sample temperature is within the boundary region for a 

oZ o Zln^ C °'T° n - ' fe Preferab,e 10 Se,6Ct b0,h 916 COrreCtk)n me,hods ad i ace m 'o the bounda j in 
order to calculate the respective correction ranges, so that a concentration of the analyte in the sample is calculated 
20 based on a value obta.nable by adding either an average value or a weighted average value to a current vat* to be 
converted into a concentration. When the environmental temperature or the sample temperature ZS^SSSi 
T^IT"' TTT* a ° CUraCy Wi " be imprWed fUrther by selectin 9 both the oo^ction methods adjacent to 

i^ss^i ssr avera9e va,ue or a wei9hted average vaiue °' 0,6 ° b,ained Je 

25 Seal iSSE""* inVemi ° n ' eXamP ' e ' 3 discri^1inan, ,unc,ion or a Mahalanobis' distance can be used as a 

rnno?! ?! fo "° Win9 explanation °o ncerns a case *> use a discriminant function as a statistical technique 
[0031] A discriminant function means a linear expression Z=f(x „ x 2> ...) defined by plural parameters (x, x, ) 
so ton™' JTT! *° ^ * * W0 P r6vious| yP r oduced groups (e.g., G1 and G2) the measurement data belong 
SL ,h.r t C ,°u!. IS < US !f ,0 discrimina,e 10 wnion of the groups the data belong. For example, a function to dis- 
criminate the level of Hot of a blood sample is defined by using a low Hot group Q, and a high Hct group G 2 as popult 

HIS- * ^ erable , t0 0btai " 3 measurement result °y P re Pann9 Plural discriminant functions in accordance with 
35 1ST a " a . £ 6 3 Sampl8, Se,eCtin9 a " 3dequate disc riminant function from the discriminant functions 
Z a nrtl T 38 a " ' ndeX ' calcu,a1ina a discriminant score Z value from the selected discriminant func- 
tion, and by correcting errors with the Z value. 

ESI .h ^ P,Ura ' dfecriminant ,unc,ions are selected by the indices, it is preferable that Z value is calculated from 
„ dominant ,u nctl0nSi and a correc , fon range ^p^g ,„ ^ z ^ „ cateulatedi so ^ me meas . 

~ t 60 fr ° m e ' ,her 3n 3Vera9e Va ' Ue ° r 3 wei9h,ed avera 9* value of the correction range. As men- 
TTJ^ *' mea * urement accuracv wi " * 'mproved by using an average value or a weighted average value of 
the correction range rather than using a correction range of any one of the functions 

ulS tha ^ reCti °K US ' n9 ^ discriminant ,unction - the c ° rrec,i °n <=an be performed to determine to which pop- 
ulation ttie data belong, by usmg the Z value code. However, it is preferable that the correction range is calculated based 

re cS ST S T ? Va ' Ue - F ° r eXample ' the discriminant score Z value can be divided into a range to be cor 
rected and a range not to be corrected (correction range is 0). 

SL .J} 6 ra " - 6 '°. be C ° rreCted 030 be divided ,Ur,her in, ° a ranae to be corrected by a correction range proper- 
oSSZl .^rri V3lUe ' and 3 range ,0 be ° 0rreC,ed bv a fixed corre o«on range without reference to the 

?nr c Z^l T aChfeVe 3 C ° rreCti0n Wftn ni9ner accuracv whe " c °mpared to a technique perform- 

so Y 3 Z™*"" 9 depend,no on ,he °' *° discriminant score. This correction technique will be regarded 

so as nonlinear correction method" in this description. ^ 

bv^electJon^ Z^Zv ° f *" T "T^ 0 "' Simi ' 3r,y> 3 measuremen ' with high accuracy can be achieved 
Lrimin^S^ h T f *?* 7* * Mana,anobis ' distance D i" P'aoe of the selection of the above-mentioned 
discnminant function and calculation of the discriminant score Z value as a statistical technique 

65 samni nh^f ""T inVen,i ° n ' e " ors b y sensor sensitivity are corrected after errors caused by 

anK E So P T?' eS aracorrec,ed - Sinc a as the parameter P2 is considered to reflect the concentration of the 

resS " ll TS^ZST T T error f 031,86(3 by me samp,e physical properties and sensor sensitivitv 

respect to the lb are corrected, and subsequently, the concentration of the analyte in the sample is calculated based on 
a value obtained by multiplying this correction value by a temperature correction coefficient 
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[0039] Furthermore, a measurement apparatus of the present invention comprises a means to measure either an 
environmental temperature or a sample temperature, a means to detect sample feeding, a means to apply a certain 
voltage at a predetermined point of time, a means to measure a current value generated by an electrochemical reaction, 
and a means to convert the measured value of the current into a concentration of an analyte in a sample, wherein the 
5 temperature measuring means, the means for detecting sample feeding, the applying means, the means for measuring 
current value and the conversion means are controlled to perform the present invention. The control is performed typi- 
cally with a microcomputer having the programmed procedures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 

[0040] 

FIG 1 illustrates a potential step in an embodiment of the present invention. FIG 2 is a flow chart to illustrate meas- 
urement and arithmetic execution in an embodiment of the present invention. FIG 3 is a detailed flow chart with 
is respect to correction I in an embodiment of the present invention. 

BEST MODE FOR CARRYING OUT THE INVENTION 

[0041 ] A typical procedure for measurement in the present invention will be explained below with a reference to an 
20 example in which the analyte is glucose in blood, and the explanation is based on FIG 1 illustrating a potential step, FIG 
2 illustrating a measurement-arithmetic execution flow, and FIG 3 illustrating a correction flow. In this example, errors in 
the terminal current value in a second excitation, which are caused by sample physical properties and sensor sensitiv- 
ity, are conected with a discriminant function, and subsequently, temperature correction is performed and the value is 
converted into a glucose concentration. 
25 [0042] FIG 1 illustrates a potential step used in the technique of the present invention. A predetermined voltage is 
applied to an electrode system of a sensor after attaching the sensor to a measurement apparatus. Current begins to 
flow when a sample is fed to the sensor. A point of time at which a current of at least a certain level is detected is deter- 
mined as a time of sample feeding, and it is defined as zero second. 

[0043] A fixed voltage is applied after the detection of the sample feeding, and the application is cancelled after t 1 
30 seconds. Therefore, an interval between 0-^ seconds is a first excitation. Either an open circuit of at most 3 seconds or 
a process to apply voltage not to cause a substantial oxidation-reduction reaction can be included just after the detec- 
tion of the sample feeding, because in the present invention, what is important is a value of maximum current or a value 
of a current thereafter in the first excitation. Voltage in the first excitation is typically in a range from 300 to 900m V, and 
tj is in a range from 3 to 15 seconds in general. As shown in the figure, the current value in the first excitation initially is 
35 raised and then lowered after having its maximum value. 

[0044] After the first excitation, either an open circuit or a voltage not to cause a substantial oxidation-reduction 
reaction is provided in the interval of t r t 2 seconds. During the interval, glucose, GOD and an electron carrier (e.g., 
potassium ferricyanide) react in a chemical reaction layer of the biosensor. The interval of t-j to t 2 is 3 to 60 seconds in 
general. During t 2 -t 3 seconds, a voltage to cause a desired oxidation-reduction reaction is applied again to the electrode 
40 in order to perform a second excitation. As a result, electrons are delivered between the electron carrier and the elec- 
trode, and a current is observed. A voltage in the second excitation is typically 300 to 900mV, and the time is typically 
2 to 1 0 seconds. The following is a series of the above-mentioned reactions. 

C 6 H l2 O 6 +GOD~»C 6 H 10 O 6 +GOD • H 2 

45 

GOD • H 2 +2Fe(CN) 6 3 -->GOD+2Fe(CN) 6 4 -+2H + 
Fe(CN) 6 4 '^Fe(CN) 6 3 -+e- 

so [0045] The above-described four parameters are calculated using a current value measured in the potential step. 
First, a current value (If (a)) of the maximum or a current thereafter in the first excitation and a current value (lb(p)) in 
the second excitation are measured to calculate P1 (lf(a)/lb(p)) and P2 (lb(p)). The If (a) is obtained typically by meas- 
uring a current value within 15 seconds after the maximum value is observed. The lb(p) preferably is obtained from a 
terminal current value 2-10 seconds after the start of the second excitation (t 2 ), and more preferably, it is obtained from 

55 the current value at the end point (t^. 

[0046] P3 (l/AI), which is a value obtained by normalizing either a derivative value or a differential value of a current 
at any point in the second excitation by the above-mentioned current value, is calculated from measured values com- 
prising a current value (lb(v» at a point (y) in the second excitation and also a current value (lb(y+A)) at a point (y+A) 
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that is positioned slightly later than the point (y) in FIG 1. The point (y) to calculate P3 is preferably in a range of 0.1 to 
1 seconds from the point of starting the second excitation (t 2 ), while A is preferably in a range of 0.1 to 1 seconds at a 
slightly later point (y+A). 

[0047] An off-value correction is performed for P3 and P4. More specifically, when the parameter value is out of an 
5 expectation range estimated from lb(P), the parameter value is substituted by a boundary value in an expectation range 
closest to the parameter value and is used for subsequent calculations. P1 and P2 are used for subsequent calculations 
without such an off-value correction. 

[0048] On the other hand, an environmental temperature or a sample temperature is measured with an internal 
thermometer in the measurement apparatus. The temperature is used in the following step to select a discriminant func- 
10 tion or to correct a calculation value by correcting errors caused by sample physical properties and sensor sensitivity. 
The environmental temperature or the sample temperature can be measured either before or after the calculation of the 
parameters. Of course, it can be measured at the same time of the calculation. 

[0049] A population is produced and used to define a discriminant function as a standard to distinguish sample 
properties and sensor sensitivity. In a measurement to produce such a population, it is essential to set conditions of a 
15 blood sample, a sensor, and a temperature. It is preferable to set as many conditions as possible for every type of sam- 
ple and concentration of the analyte. At least two types of groups, e.g., a high value and a low value are set with respect 
to the sample physical properties. Specifically, measurement data based on a high Hct and a low Hct will be taken when 
the condition is Hct. 

[0050] At least 10 lots of sensors are selected so that the manufacture season (elapsed time after manufacture), 
20 the manufacture line, sensor sensitivity or the like correspond to a general state of use as much as possible. As many 
temperature settings as is practical are preferred, and at least two types of settings are required with a boundary of 30 
to 35°C, since reactivity of If/lb used as the parameter with regard to the concentration of the analyte is shaped like a 
quadratic function having a downward convex curve if the temperature is equal to or lower than the above-mentioned 
boundary, while the reactivity will be shaped like a substantially straight line if the temperature is equal to or higher than 
25 the boundary. 

[0051] The population is produced by combining the sample and the sensor, which are measured under the respec- 
tive temperature conditions. Namely, a correction method to deal with errors caused by the sample properties and the 
sensor sensitivity is prepared under the set conditional temperature. In a correction method explained later in this 
description, two types of corrections, that is, a lower-and-medium-temperature side correction and a high-temperature 

30 side correction are included, in which the boundary is 30 to 35°C. 

[0052] First, the terminal current value lb(p) in the second excitation is corrected in the lower-and-medium-temper- 
ature side correction. As shown in FIG 3, this correction includes a correction of errors caused by the sample properties 
and a correction of errors caused by the sensor sensitivity, which are performed in this order. 
[0053] Discriminant functions for correcting sample physical properties are prepared for the respective ranges of 

35 the analytes' concentrations. The discriminant functions for correcting sample physical properties are produced as dis- 
criminant functions to use the parameters for discriminating to which of two populations a certain measurement data 
belongs, where the populations are classified depending on the heights of Hct. Preferably in this case, data including a 
parameter considered as an off-value is excluded from the populations, since the populations as standards are pre- 
ferred to be uniform. 

40 [0054] The concentration range of the analytes can be set experientially. 

[0055] The discriminant function(s) for correcting sample physical properties can be selected according to an index 
of a value multiplied by lb(P); it/they can be selected according to an index of a value obtained from a temperature cor- 
rection coefficient (T1) multiplied by lb(P), where the temperature correction coefficient is previously calculated from an 
environmental temperature or a sample temperature measured with a thermometer. 

45 [0056] The following is an explanation about a case where a single discriminant function for correcting sample 
physical properties is selected. First, a parameter calculated in the above-mentioned manner is substituted for the dis- 
criminant function for correcting sample physical properties so that a discriminant score Z value is calculated. It is pos- 
sible to determine to which population the data belongs to, based on this Z value code. Or the lb(P) can be corrected 
by determining a correction range from the relation between lb(P) and the Z value. For example, the Z value is classified 

so into five regions as described below, and the correction range is determined with a nonlinear correction. In the following, 
"correcting with a fixed range" means correcting with a fixed range determined by a measured lb(P), regardless of the 
Z value. The fixed range can be determined experientially by taking the types of the analytes or some other factors into 
consideration. 

[0057] A linear correction means to determine a correction range from a proportional relation with the Z value with 
55 an upper limit of the above-mentioned fixed range. 

[0058] Thus determined correction range is regarded as a correction range A. The five regions in the example can 
be determined experientially taking the types of the analytes and some other factors into consideration. 
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(1) When the Z value is positioned toward the extreme plus side, it is corrected toward the minus side with a fixed 
range. 

(2) When the Z value is positioned toward the plus side, it is linear-corrected toward the minus side. 

(3) When the Z value is positioned in the vicinity of zero, no correction is performed. 

5 (4) When the Z value is positioned toward the minus side, it is linear-corrected toward the plus side. 

(5) When the Z value is positioned toward the extreme minus side, it is corrected toward the plus side with a fixed 
range. 

[0059] When errors caused by sample physical properties is not corrected (correction range A=0), errors caused 
w by sensor sensitivity are corrected. Discriminant functions for correcting sensor sensitivity also are prepared for the 
respective concentration ranges of anafyte. The discriminant functions for correcting sensor sensitivity are selected with 
an index of a value obtained by multiplying by lb(p). When a single discriminant function for correcting sensor sensitivity 
is selected, a discriminant score Z value is calculated from the parameter and the discriminant function, and a correc- 
tion range corresponding thereto is calculated. The discriminant function for correcting sensor sensitivity is produced 
is by using the parameter, as a discriminant function to discriminate populations classified based on differences in the 
electrode sensitivity or the like. When plural discriminant functions for correcting sensor sensitivity are selected, discri- 
minant score Z values and correction ranges corresponding to the Z values are calculated from the respective discrimi- 
. nant functions for correcting sensor sensitivity, and the average value of the correction ranges is determined to be a 
correction range B. 

20 [0060] When the above-mentioned discriminant function(s) for correcting sensor sensitivity is/are used, a correction 
range intermediate value C is determined in the following manner. 

(1) When correction range A*0, 

correction range intermediate value C = correction range A 
25 (2) When correction range A=0, 

correction range intermediate value C = correction range B 

[0061 ] When plural discriminant functions for correcting sample physical properties are selected, a correction range 
intermediate value C is calculated in the above-mentioned manner from the respective discriminant functions for cor- 

30 recting sample physical properties, and the average value or a weighted average value is determined as a final correc- 
tion range D. When a single discriminant function for correcting sample physical properties is selected, the correction 
range intermediate value C is the final correction range D. A value obtained by adding the final correction range D to 
the terminal current value (lb(P» in the second excitation is the final current value lb(p) lower-middle )in the lower-and- 
middle-temperature side correction. 

35 [0062] In a high-temperature side correction, procedures such as data-taking, production of populations and discri- 
minant functions and calculation of correction range are performed similar to the lower-and-middle-temperature side 
correction, so that a high-temperature side final current value lb(P) high>is calculated. 

[0063] From the temperatures measured with a thermometer, a pre-temperature correction current value !b(P)' is 
obtained. For example, lb(P)'=lb(P) flower-middle > at a temperature from 10 to 30°C, and lb(P)'=lb(P) high ) at a tem- 

40 perature from 35 to 40°C. The temperature range from 30 to 35°C is determined as a boundary region, and either the 
average value or the weight average value of lb(p) tower-middle >and lb(p)high)is determined as lb(P)'. The boundary 
region can be selected experientially by taking the types of the analytes and some other factors into consideration. 
[0064] The lb(p)' is multiplied by a temperature correction coefficient (T1 ) that is calculated previously from an envi- 
ronmental temperature or a sample temperature. The thus obtained value is determined as a post-correction final cur- 

45 rent value lb(p)" of lb(P). And the lb(P)" is converted into a glucose concentration, using a calibration curve or a 
calibration formula that is produced previously. 

[0065] In the respective sections to select discriminant functions for correction as shown in FIG. 3, only ® (Q)' or 
© M for a discriminant function for correcting sensor sensitivity) is selected if a single type should be selected. The cor- 
rection ranges are calculated respectively from ® and © ((J)' and or © H and ©" for discriminant functions for 
so correcting sensor sensitivity) if two types are selected, and the average values are determined as the correction range 
intermediate value and the final correction range. 

[0066] The present invention will be explained in detail with respect to specific embodiments. 
[0067] A glucose sensor system to quantify a glucose concentration in a blood sample will be taken below as an 
example of a biosensor system. 
55 [0068] A glucose sensor used in the glucose sensor system is configured as follows. On an insulating substrate 
comprising polyethylene terephthalate (PET), a carbon electrode system and an electrical insulating layer are formed 
by screen printing, and the carbon electrode system comprises a silver lead, a working electrode and a counter elec- 
trode. The electrical insulating layer has a predetermined area in which the working electrode and the counter electrode 
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are exposed, and partially covers the silver lead. On the thus formed electrode system, a layer of carboxymethylcellu- 
lose (CMC) as a hydrophilic polymer and further, an (enzyme+mediator+CMC) layer are formed. The (enzyme+media- 
tor+CMC) layer comprises glucose oxidase (GOD) as an enzyme, potassium ferricyanide as a mediator (electron 
carrier), and a CMC. The combination of the CMC layer and the (enzyme+mediator+CMC) layer is regarded as a chem- 

5 ical reaction layer. An insert comprising a cover and a spacer further is formed, so that a fixed amount (about 3\lL) of a 
blood sample contacting with the insert is fed by a capillary phenomenon to the chemical reaction layer and to the elec- 
trode system. For the measurement apparatus, a commonly-used potentiostat 100B/W (supplied by BAS) was used. 
While a typical biosensor system uses a combination of a single test cell and an exclusive apparatus, a generalized 
apparatus is used in the present embodiments. However, the embodiments can be performed easily even if a simple 

to apparatus of a conventional technique is used. The measurement in the embodiments was performed in an environ- 
mental test lab, but the embodiments can be performed also by using a simple apparatus containing a temperature sen- 
sor. 

[0069] The potential steps were provided in the following manner. 

[0070] A glucose sensor was attached to an apparatus, and a voltage of 500mV was applied to the glucose sensor 
is before sample feeding and for 7 seconds after the detection of sample feeding. After an open circuit of 18 seconds, a 
voltage of 500mV was applied again for 5 seconds. The voltage was applied each time for the following purposes. The 
application before the sample feeding is to detect a current generated electrochemically by the sample feeding. Sample 
feeding is detected as a current of at least a certain level, and thus, t=0 will be defined. The 7-second application after 
the sample feeding is for the first excitation, and the main purpose is a detection of difference in the current profile 
20 caused by the difference in the sample physical properties. The open circuit of 18 seconds is to promote the chemical 
reaction and to store Fe(CN) 6 4 * as a matrix for an electrochemical reaction. The subsequent 5-second application of 
voltage is for the second excitation, and the main purpose is to convert the Fe(CN)6 4 " to Fe(CN)6 3 *electrochemically, 
and to detect a difference in the current profile caused by the difference in the sensor sensitivity and the concentration 
of the analyte's concentration. 

25 

(First Embodiment) 

[0071] First, errors caused by Hct were corrected, since sample physical properties are the most influential factor 
to cause errors in a measurement result. 

30 [0072] The current was measured at every 0. 1 seconds just after the detection of the sample feeding and also after 
the re-application, and current values at the following three points were used for calculation: the peak value l p that was 
measured just after the detection of the sample feeding; 1(0.5) that was measured 0.5 seconds after the re-application 
of current; and 1(5) that was measured 5 seconds after the re-application of current The l p was calculated in the follow- 
ing manner. Namely, l p was l(t) at the moment of l(t)>l(t+0.1) in a current just after the detection of the sample feeding. 

35 The respective current values were measured for the following purposes. The l p was measured to detect a difference 
in the current profile caused by a difference in the sample physical properties such as a blood viscosity. The 1(0.5) was 
measured to see the ratio to l p , so that the glucose concentration factor of l p would be reduced. The 1(5) was measured 
to obtain a current value that depends on the glucose concentration and that can be used for a glucose concentration 
conversion. 

40 [0073] l(0.5)/l p was obtained for an index to indicate the blood viscosity, and it was determined as ce^ In this embod- 
iment, the following experiment was performed indicating erythrocyte amount with hematocrit, since erythrocyte is a . 
blood cell in a blood sample and it can be adjusted easily. 

[0074] A relation between a t and 1(5) is shown below, and the relation distinguishes a specimen of 30% hematocrit 
and a specimen of 60% hematocrit regardless of the glucose concentration in the blood sample. 

45 ' 

z=-23.1 15xa 1 +1 .254xl(5)+41 .580 

[0075] The level of hematocrit in the blood sample is reflected by the value of z. When z is positively large, the 
hematocrit in the blood sample is considered to be small. 
so [0076] A current value l'(5) after the correction of the sample physical properties was calculated by the following 
expression from the relation between a 1 and 1(5). 

l'(5)=slopex(-z)+l(5) 

55 In this equation, slope=0.0238xl(5)-0.0338. 

[0077] The thus calculated current value l'{5) with a corrected blood viscosity will be converted into a glucose con- 
centration value by using a certain calibration curve table. 

[0078] The following data were obtained using a blood sample (glucose viscosity: 536mg/dL, hematocrit: 30%). 
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l p =10.210(uA) 
l(0.5)=20.290(nA) 

5 l(5)=10.110(|iA) 

[0079] Regarding the data, a current value l*(5) with a corrected blood viscosity can be obtained through the follow- 
ing arithmetic executions. 

w a 1 =l(0.5)/l p =1.987 

z=-23.1 15xa ,+1 .254xl(5)+41 .58=8.328 

slope=0.0238xl(5)-0.0338=0.207 

15 

!'(5)=slopex(-2)+l(5)=8.387 

This l'(5) will be converted into a glucose concentration value by using a calibration curve that has been produced. 
[0080] A blood sample used in this embodiment was prepared in the following manner. Blood samples collected 
20 from four persons with the identical blood type were mixed to prepare blood samples having hematocrit of 30%, 45% 
and 60%, respectively. An aqueous solution of glucose was added to these blood samples in order to adjust the glucose 
concentrations to 83mg/dL, 307mg/dL, and 536mg/dL. Measurement was performed using these 9 types of blood sam- 
ples, where the measurement number n was determined to be 15. 

[0081] As explained later in detail, errors caused by hematocrit with regard to blood samples were corrected suc- 
25 cessfully in this embodiment. 

[0082] A coefficient of variation (C.V. (%)) is obtained by a calculation based on a group of measurement data for 
blood samples having different hematocrits for the respective glucose concentrations. The following Table 1 shows the 
coefficient of variation to compare the values before and after the correction of the sample physical properties. 

30 

Table 1 


Glucose concentration 

Pre-correction 

Post-correction 

83mg/dL 

5.40% 

3.92% 

307mg/dL 

10.69% 

5.67% 

536mg/dL 

16.52% 

8.56% 


[0083] The hematocrit differs depending on individuals. When the standard value is 45% in the present invention, 
40 the divergence (%) of the hematocrit 45% specimens from the measured values for the respective glucose concentra- 
tions are as shown in Tables 2 and 3. Here, divergence of the post-correction data means a divergence from a corrected 
average value of hematocrit 45%. Table 2 shows the divergences (%) from measured average values with respect to 
the cases of hematocrit 30% and of hematocrit 45%. 

45 

Table 2 


Glucose concentration 

Pre-correction 

Post-correction 

83mg/dL 

4.27% 

-0.15% 

307mg/dL 

14.43% 

1.50% 

536mg/dL 

15.69% 

-1.92% 


[0084] It is clear that errors caused by the physical properties of the blood samples were improved as a result of the 
55 correction. 

[0085] Table 3 shows divergence (%) of measured average values for hematocrit 60% from those for hematocrit 
45%. 
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Table 3 


Glucose concentration 

Pre-correction 

Post-correction 

83mg/dL 

-4.31% 

-1.77% 

307mg/dL 

-10.85% 

-4.10% 

536mg/dL 

-22.21% 

-9.12% 


10 

[0086] It is clear that errors caused by the physical properties of the blood samples were improved as a result of the 
correction. 

(Second Embodiment) 

15 

[0087] In this embodiment, a correction of errors caused by an influence of Hct as sample physical properties and 
a correction of errors caused by sensor sensitivity are performed simultaneously in a series of arithmetic executions. 
[0088] The current was measured for every 0.1 seconds after the detection of sample feeding, and the following five 
points were used for the executions: a value If (7) that was measured 7 seconds alter the first excitation, a value lb(0.3) 
20 that was measured 0.3 seconds after the second excitation, a value lb(0.5) that was measured 0.5 seconds after the 
second excitation, a value lb(1) that was measured 1 second after the second excitation, and a value lb(5) that was 
measured 5 seconds after the second excitation. 
[0089] The following four parameters were used for the operation: 

25 P1:lf(7)/lb(5) 
P2:lb(5) 

P3: l/AI(0.5)=lb(0.5)/{lb(0.5)-lb(1 )} 
P4:lb(0.3)/Jb(5). 

30 [0090] Regarding l/Al(0.5) and lb(0.5)/lb(5), an off-value correction was performed for the data deviating from a 
range expected from !b(5) that is considered as an index for a glucose concentration. Specifically, when the parameter 
value was higher than the upper limit of the range expected from Ib(5), the upper limit was adapted as the parameter. 
When the value was lower than the lower limit of the range expected from lb(5), the lower limit was adapted as the 
parameter. A range expected from lb(5) was adapted as a median, and the parameter range was set as the median 

35 ±2SD from the data derived from a population. For example, when the range of l/AI(0.5) expected from the measured 
lb(5) was 3.68 and the standard deviation of the parameter was 0.39, a range expected from the lb(5) was from 2.90 to 
4.46, and if the actually obtained l/AI(0.5) was 4.50, it was considered as l/AI{0.5)=4.46 . 

[0091 ] Data for making populations for discriminant analysis were obtained at environmental temperatures of 25°C 
and 32°C. The blood samples were prepared by mixing blood collected from plural subjects and by adjusting the glu- 
40 cose concentrations to be 30, 60, 75, 100, 150, 200, 250, 300, 400, and 500mg/dL, and Hct to be 25, 45, and 65%. A 
glucose sensor manufactured within a year and a half was used, and 10 to 14 lots were selected for every sample so 
that the lots are selected at random regardless of the manufacture season, manufacture line, and sensor sensitivity. In 
this embodiment, the following two populations were prepared for every environmental temperature: (1) a population for 
a Hct correction as a correction of sample physical properties, and (2) a population for correcting sensor sensitivity. 
45 Regarding (1), a group of integrated data of three types of Hct specimens with the respective glucose concentrations 
were classified evenly into five groups depending on the value of lb(5), and a group (H,) with the largest lb(5) and a 
group (L^ with the smallest lb(5) were adapted as populations. Populations were prepared in this manner without sim- 
ply taking groups with Hct 25% and with Hct 65%, since correction of the sensor sensitivity was performed preferably 
as well as the correction of Hct. Regarding (2), a group of data specimens with Hct 45% for the respective glucose con- 
so centrations were classified evenly into'five groups depending on the value of lb(5), and a group (H 2 ) with the largest 
lb(5) and a group (L 2 ) with the smallest lb(5) were adapted as populations. Here only the data group with Hct 45% was 
adapted to improve accuracy in correcting sensor sensitivity, since 45% was considered as a standard Hct value. 
[0092] A discriminant function is produced as a function to discriminate the populations that have been obtained 
using parameters for the respective measurement temperatures. Specifically, a discriminant function for Hct correction 
55 discriminates between H, and , and a discriminant function for correcting sensor sensitivity discriminates between H 2 
and Lg. In the following, correction I is a correction method using an discriminant function obtained from a measurement 
data at 25°C, while correction II is a correction method using a discriminant function obtained from a measurement data 
at 32°C. The correction I indicates a correction under a condition of a lower or room temperature, while the correction 


10 


EP 1 081 490 A1 


li indicates a correction under a high temperature condition. The main reason for classifying the correction method into 
two categories depending on the temperature is that If (7)/lb(5) is considerably temperature dependent. 
[0093] Four types of discriminant functions for Hct correction with regard to the correction I were produced: func- 
tions for 30-75mg/dL, 75-1 50mg/dL, 1 50-300mg/dL t and 300-500mg/dL In a case of 75-1 50mg/dL, a discriminant func- 

5 tion to discriminate two populations of Hj and L t for 75mg/dL and 1 50mg/dL was obtained. In this case, the discriminant 
function was produced so that the discriminant score Z value of a data having high values due to the influence of Hct 
(i.e., a data of high Hct specimens) would be plus. In producing the populations, data including parameters regarded as 
abnormal values as a result of a null hypothesis test of (Grubs-Smirnor were eliminated, since the populations were 
preferably as uniform as possible. 

w [0094] The discriminant function for Hct correction was produced by selecting a value as an index obtained by mul- 
tiplying a temperature correction coefficient (T1) by lb(5), where the temperature correction coefficient was calculated 
previously from the environmental temperature, and the selection numbers were 1-2. f 1 was determined as 1 at 25°C, 
a temperature that is considered as a standard. A case where plural discriminant functions were selected will be 
described later. When a single discriminant function for Hct correction was selected, the parameter was substituted for 

15 the discriminant function, and a discriminant score Z value and a correction range A corresponding thereto were calcu- 
lated. The correction range A was determined with a nonlinear correction method. Namely, when 
Z A <Z B <Z C <Z D (Z A ,<Z B ,<Z C ,<Z D are constants experientially determined for every discriminant function) and when a 
current value lb(5) and a discriminant score Z value were obtained, corrections were performed as follows: when Z<Z A , 
correction was performed with a fixed range toward the plus side regardless of Z value; when Z A <Z<Z B , a linear correc- 

20 tion was performed toward the plus side corresponding to the Z value; when Z B <Z<Z C , no correction was performed; 
when-Z c <Z<Zr> a linear correction was performed toward the minus side corresponding to the Z value; and when Zq<Z, 
a correction with a fixed range was performed toward the minus side regardless of the Z value. The correction with a 
fixed range means correcting with a fixed range determined by a measured lb(5) for every discriminant function. A linear 
correction is determining a correction range from a proportional relation with the Z value with an upper limit of the 

25 above-mentioned fixed range. For example, if a linear correction is performed toward the minus side when 5<Z<1 0, the 
correction range is -0.10 for Z>10, a correction range for Z=8 is calculated according to the following equation (1). 

Correction range A=-0.10x(8-5)/(10-5)=-0.06 (1) 

30 [0095] The term 'correction of a fixed range' means that a correction range will have a fixed value once lb(5) is 
determined. The correction range itself varies in accordance with lb(5). When no Hct correction was performed (correc- 
tion range A=0), the sensor sensitivity as mentioned above was corrected as described below in order to obtain a cor- 
rection range B. When Hct correction was performed (correction range A*0), correction of sensor sensitjvity will not be 
performed. 

35 [0096] Seven types of discriminant functions for correcting sensor sensitivity with regard to correction I were pre- 
pared: functions for 75-100mg/dL, 100-150mg/dL, 150-200mg/dL, 200-250mg/dL, 250-300mg/dL, 300-400mg/dL and 
400-500mg/dL In a case of 75-100mg/dL, a discriminant function to discriminate two populations of H 2 and L 2 for 
75mg/dl_ and 1 0Omg/dL was obtained and the function was determined as a discriminant function for correcting sensor 
sensitivity for the above-identified concentration ranges. In this case, the discriminant function was produced, so that 

40 the discriminant score Z value measured with a highly sensitive sensor would be plus. In producing the populations, 
data including parameters regarded as extremely high/low as a result of a null hypothesis test of (Grubbs-Smirnov were 
eliminated, since the populations were preferably as uniform as possible similar to the above-mentioned case of pro- 
ducing a discriminant function for Hct correction. 

[0097] The discriminant function for correcting sensor sensitivity was selected with lb(5) as an index, and the selec- 
45 tion numbers were 1-2. When a single discriminant function was selected, the parameter value was substituted for the 
discriminant function, and a discriminant score Z value and a correction range corresponding thereto were calculated. 
The correction range was regarded as a correction range B. The correction range was determined with a nonlinear cor- 
rection method similar to the case of Hct correction. When plural discriminant functions were selected, correction 
ranges were calculated from the respective discriminant functions, and the average value was regarded as the correc- 
so tion range B. 

[0098] A correction range intermediate value C was determined as follows when a single discriminant function for 
Hct correction was selected. 

(1) When correction range A*0, correction range intermediate value C=correction range A 
55 (2) When correction range A=0, correction range intermediate value C=correction range B 

[0099] When two types of discriminant functions for Hct correction were selected, a correction I final correction 
range D was obtained in a similar way by obtaining an average value or a weighted average value of the correction 
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range intermediate values C obtainable from the respective discriminant functions for Hct correction. When only one 
type of discriminant function for Hct correction was selected, a correction range intermediate value C thereof was deter- 
mined as the correction I final correction range D. The correction I final current value lb(5)| was a value obtained by add- 
ing the final correction range to lb(5). 
5 [01 00] Regarding the correction II, the technique of the correction I was adapted to take data, produce populations 
and discriminant functions, and to calculate correction ranges, so that a Correction II final current value lb(5)n was . 
obtained. 

[0101] A pre-temperature correction current value !b(5)* was determined as follows from a temperature T (°C) 
detected with a temperature sensor. 

10 

WhenT<30, lb(5)'=lb(5) , . 

When SOiTT^S, lb(5)'=(35-T)/(35-30)xlb(5) ,+(T-30)/(35-50)xlb(5) „ . 
WhenT>35, lb(5)'=lb(5)„. 

is [01 02] For example, lb(5)' can be obtained from the following equation when T=32°C: 

lb(5)*=3£xlb(5) ,+2/5xlb(5) ,, 

[0103] lb(5) n was obtained by multiplying thus obtained 1b(5)' by a temperature correction coefficient (T2). This 
20 lb(5) M was converted into a glucose concentration value according to a calibration curve table that was produced sep- 
arately. 

[01 04] For example, the following data were obtained by measuring blood specimens of a glucose concentration of 
154mg/dL and Hct 65% under an environmental temperature of 25°C. 

25 lf(7)=2.142(nA) 

lb(0.3)=7.578{uA) 

lb(0.5)=6.312(nA) 

30 

lb(1)=4.932(nA) 
lb(5)=2.142(uA) 

35 [01 05] The following parameters are obtained from the above-identified data. 

If(7)/lb(5)=0.8069 
l/Al(0.5)=lb(0.5)+{lb(0.5)-lb(1)}=4.575 

40 

lb(0.3)/lb(5)=2.855 
lb(5)=2.654 

45 [0106] At this time, expectation ranges of l/AI(0.5) and lb(0.3)/lb(5) are calculated from the value of lb(5) respec- 
tively to be 3.050-4.496 and 2.817-4.361. Since l/AI(0.5) exceeds the upper limit of the expectation range and thus is 
considered as an off-value, l/AI(0.5)=4.496 in the subsequent calculations. 

[01 07] Correction I is performed. A discriminant function for Hct correction is selected with an index of lb(5)=2.654. 
When lb(5) ranges from 2.407 to 3.400, a discriminant function for 75-150mg/dL and a discriminant function for. 150- 
50 300mg/dL are used. The discriminant functions are defined respectively as follows. A discriminant function for 75- 
1 50mg/dL • Hct correction: 

Z(L)=34.525xP1-27.333xP2-21.717xP3+3.1870xP4+141.90. 

55 A discriminant function for 150-30mg/dL* Hct correction: 

Z(N)=331 .25xP1 -1 6.739xP2-22.331 xP3+6.71 91 xP4-1 58.00. 
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In these equations, P1 =lf{7)/lb(5) , P2=I/AI(0.5) . P3=lb(0.3)/lb(5) , and P4=lb(5) . 

[0108] In this embodiment, Z(L)=-5.753. Since a linear correction is performed toward the plus side when - 
8<Z(L)<-4, a correction range is calculated from the following equation: 

A L ={0.1 205xlb(5)-0. 1 049}x{-4-Z}/{-4-(-8)}=0. 1 024. 

[0109] The sensor sensitivity will not be corrected since A^O. The value of 0.1024 will be a correction range inter- 
mediate value C L . 

[01 10] . In this embodiment, Z(N)=-1 1 .887 . Since no correction is performed when -1 5<Z(N)<5, the sensor sensitiv- 
ity will be corrected. Since lb(5)xtemperature correction coefficient (T1 )=2.654 and this figure is within a range of 2.582- 
3.316, a discriminant function used for correcting sensor sensitivity is of 150-200mg/dL. The discriminant function of 
150-200mg/dL is defined as follows. A discriminant function for correcting sensor sensitivity of 1 50-200mg/dL: 

Z(S 3 )=43.581 xP1 -25.041 xP2-25.472xP3+2.0536xP4+1 41 .86. 

In this equation, P1=lf(7)/ib(5) , P2=I/AI(0.5) , P3=lb(0.3)/lb(5) , and P4=lb(5) . 
[01 1 1 ] In this example, Z(S 3 )=-2.827 . 

[0112] Since a linear correction is performed toward the plus side when -6<Z(S^<-2, a correction range B N is cal- 
culated from the following equation. This B N will be a correction range intermediate value C N . 

B N =C N ={0.0312xlb{5)-0.0121}x{-2-Z(S3)}/{-2-(-6)}=0.0146 

[0113] Since two discriminant functions for Hct correction are selected in this example, a final correction range D, 
of the correction I is an average value obtained from the two correction range intermediate values C L and C N : 

D ,=(0.1 024+0.01 46)+2=0.0585. 

Therefore, a post-correction current value lb(5) t obtainable by the correction I can be obtained in the following manner. 

Ib(5) ,=2.654+0.0585-2.713 

[0114] Since the temperature is 25°C and the pie-temperature correction current value lb(5)'=lb(5) , , calculation 
procedures for obtaining lb(5)n will be omitted from the following description. Since the temperature is 25°C, the tem- 
perature correction coefficient is 1 and the post-temperature correction current value lb(5) ,, =lb(5)'=2.713 . This lb(5)" 
value is substituted for a calibration curve produced separately in order to be converted into a glucose concentration. 
An average value of lb(5) of a data of a glucose concentration of 154mg/dL and Hct 45%, i.e., an expectation value of 
lb(5) for a glucose of 154mg/dL was 2.790. A correction according to the present invention decreased divergence from 
the expectation range (pre-correction value=2.142; post-correction value=2.713). 

[0115] In this embodiment, errors caused by sample physical properties and sensor sensitivity were decreased. 
This will be explained below referring to an result of a case where the glucose concentration was 150mg/dL 
[0116] Ten lotsxnumber N (10) of sensors were measured for every Hct specimen. The coefficient of variation 
(CV)(%) for the entire measurement data varied before and after the correction as described below. The table below 
shows that errors caused by sensor sensitivity is decreased. 


Table 7 


Specimen Hct 

Pre-correction 

Post-correction 

Correction effect (Pre- 
correction- Post-correc- 
tion) 

25% 

4.7% 

4.3% 

A0.4% 

45% 

4.5% 

4.2% 

A0.3% 

65% 

5.1% 

4.7% 

AO.4% 


[0117] Since the hematocrit differs depending on individuals, the standard value is 45% in the present invention. 
Divergences (%) of the respective Hct specimen current average values from a Hct 45% specimen current average 
value are shown in Table 2. It is clear from this table that errors caused by Hct were corrected further. 
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Table 8 


Specimen Hct 

P re-correction 

Post-correction 

Correction effect (Pre- 
correction-Post-correc- 
tion) 

25% 

15.6% . 

8.2% 

A7.4% 

65% 

-7.2% 

-3.8% 

A3.4% 


[01 18] As an index to indicate a maximum divergence of an obtained measurement value from a true value or from 
an expectation value, a total error (TE) is considered. TE is defined as follows. 

15 TE(%)=DEV(%)+1.65xCV(%) 

[01 19] A coefficient 1 .65 of a CV term is a value of Z value (a value not as a discriminant score but as a value in a 
standard normal distribution) composing one-sided probability of 5% of a standard normal distribution, and TE defined 
in this term is a divergence range from a true value to which 95% of the entire data is considered to be included. 
20 [0120] TE in this example is as shown in Table 3. 


Table 9 


25 


30 


Specimen Hct 

Pre-correction 

Post-correction 

Correction effect (Pre- 
correction-Post-correc- 
tion) 

25% 

23.3% 

15.3% 

A8.0% 

45% 

7.4% 

7.0% 

AO.3% 

65% 

15.6% 

11.5% 

A4.4% 


INDUSTRIAL APPLICABILITY 

35 [01 21 ] As mentioned above, according to the measurement method of the present invention, errors with respect to 
analytes caused by factors such as sample physical properties and sensor sensitivity can be corrected efficiently with- 
out modifying compositions of a chemical reaction layer or an electrode system of a sensor, and thus, a reliable meas- 
urement value can be obtained. 

40 Claims 

1 . An electrochemical measurement method to measure a concentration of an analyte in a sample by using a biosen- 
sor, comprising calculation of parameters from a current value that is obtained as a result of application of a fixed 
voltage to the biosensor after sample feeding, and correction of errors of a measured value using these parameters 

45 by a statistical technique. 

2. The measurement method according to claim 1 , wherein parameters that are affected by errors caused by sample 
physical properties are selected. 

so 3. The measurement method according to claim 1 , wherein parameters that are affected by errors caused by sensor 
sensitivity are selected. 

4. The measurement method according to claim 1 , wherein the concentration of the analyte is obtained by applying a 
predetermined voltage twice to the biosensor after sample feeding, calculating parameters P1 and P2 from the cur- 
55 rent values provided by the application, and by correcting the errors from these parameters by a statistical tech- 
nique, where 

P1 : a ratio (If/lb), where (If) is the value of maximum current or a current occurring after the maximum in a first 
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excitation, and (lb) is the value of current read at any point in a second excitation; 
P2: a current (lb) read at any point in the second excitation. 

5. The measurement method according to claim 4, wherein the current value lb at any point in the second excitation 
5 used for the parameters P1 and P2 is a terminal current value (lb(P)) in the second excitation current. 

6. The measurement method according to claim 4, wherein at least one of parameters P3 and P4 is used together 
with the parameters P1 and P2, where 

10 P3: a value (l/AI(5)) obtained by normalizing a derivative value or a differential value at any point in the second 

excitation with a current value at the same point; and 

P4: a ratio (lb(a)/lb(P)) of an initial current value (lb(ct)) to a terminal current value (lb(P)) in a second excitation. 

7. The measurement method according to claim 6, wherein at least one of the parameters P3 and P4 is substituted 
15 with a closest expectation value when the parameter is out of an expectation range. 

8. The measurement method according to claim 1, wherein an adequate correction is performed in accordance with 
. an environmental temperature or a sample temperature by preparing a plurality of correction formulas using a sta- 
tistical technique corresponding to the environmental temperature or the sample temperature and by selecting an 

20 optimum correction formula from the plural correction formulas. 

9. The measurement method according to claim 8, wherein the concentration of the analyte in the sample is obtained 
when the environmental temperature or the sample temperature corresponds to the boundary region of the tem- 
perature range used for the correction formula, by selecting both the correction formulas adjacent to the boundary, 

25 calculating correction ranges from the respective formulas, adding the average values or weighted average values 
to the current value for concentration conversion, and obtaining the concentration depending on the additional 
value. 

10. The measurement method according to claim 1, wherein a discriminant function is selected from the parameters. 

30 

11 '. The measurement method according to claim 10, wherein a discriminant score is obtained from the discriminant 
function and the concentration of the analyte is calculated from a correction range proportional to the discriminant 
score. 

35 12. The measurement method according to claim 10, wherein a range to calculate the concentration of the analyte 
depending on a fixed correction range regardless of the discriminant score, and a range to calculate the concen- 
tration of the analyte without performing a correction are set. 

1 3. A measurement apparatus comprising a meter for measuring an environmental temperature or a sample tempera- 
40 ture, a detector for detecting sample feeding, a device for applying a certain voltage to a biosensor at a predeter- 
mined point of time, a meter for measuring a value of current generated by an electrochemical reaction, and a 
converter for converting the measured value of the current into a concentration of the analyte in the sample; 
wherein the meter for measuring temperature, the detector for detecting sample feeding, the device for application, 
the meter for measuring the current value and the converter are controlled to perform the measurement according 
45 to any of claims 1 to 1 2. 
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